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^ 1. Motivation 

^ A limited number of instabilities are responsible for most of 
5^, the pattern formation in the Universe. Stars are formed by the 
Jeans instability, and the heat transport within them is driven 
by the Rayleigh-Benard instability. Most high-energy radia- 
tion, however, is produced by disks around compact objects 
and black holes. These accretion disks must therefore be tur- 
bulent. It is now believed that this turbulence is caused by 
the instability of their Keplerian rotation law in the presence 
of weak large-scale magnetic fields. This 'magnetorotational 
instability' (MRI) was first discovered by Velikhov (1959), 
who considered an electrically conducting fluid between ro- 
tating cylinders with J7out < ^in- The existence of this insta- 
bility has been confirmed analytically and numerically many 
times; see for example Balbus (2003) and Riidiger & Holler- 
bach (2004) for recent reviews. 

There is considerable interest in studying the magnetoro- 
tational instability in the laboratory. The simplest, and most 
widely studied design, is to confine a liquid metal between 
differentially rotating cylinders, and impose a magnetic field 
along the axis of the cylinders. This yields the MRI for mag- 
netic Reynolds numbers of order 10. However, because of 
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the extremely small magnetic Prandtl numbers Pm(= i^/i]) 
of available liquid metals (see Table QJ, that translates into 
ordinary Reynolds numbers exceeding 10^, which causes se- 
vere difficulties (HoUerbach & Fournier 2004). In order to 
overcome such problems, HoUerbach & Riidiger (2005) pro- 
posed imposing an azimuthal field as well, with dramatic con- 
sequences, namely a reduction of the critical Reynolds num- 
ber Rec from O(IO^) to 0(10'^). These new solutions are 
also essentially independent of Pm, and are therefore ideally 
suited to experimental realisations in the laboratory. These re- 
sults were for insulating inner and outer cylinders, in which 
case the axial current one must impose to generate this az- 
imuthal magnetic field is around 2500 A (using liquid sodium 
as the fluid), which is close to the upper limit of what is ex- 
perimentally possible. Here we therefore consider conducting 
cylinders, and show that the necessary current (as well as the 
Reynolds number) is smaller than for insulating cylinders. 



Table 1. Liquid metal material parameters in c.g.s. 
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J? 


Pm 




sodium 
gallium 


0.9 
6.0 


7.1 
3.2 


IQ-^ 
IQ-^ 


0.8 • 10^ 
2.1 • 10^ 


0.9- 10"^ 
1.5 • lO"*^ 


8.15 
22.0 



3 THE EQUATIONS 



2. Introduction 

We consider the stability of the flow between two rotating 
coaxial infinitely long cylinders, in the presence of a constant 
axial field B^, and a current-free (within the fluid) azimuthal 
field Bij). Figure n shows a sketch of the geometry. The fluid 
is incompressible, with density p, kinematic viscosity v and 
magnetic diffusivity rj. 




Fig. 1. (online colour at www.an-journal.org) The basic ge- 
ometry of the problem, consisting of two cylinders of radii 
Rin and i?out, rotating at and l?out- Bz and B^ are the 
externally imposed magnetic fields, B^ by an axial current 
inside the inner cylinder. 



By conservation of angular momentum the basic state dif- 
ferential rotation profile is given by 



r2 = a+^, 

where a and b are given by 

^2 1 
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with fi = r2out/An and r) = i?in/i?out- -Rin and i?out are the 
radii of the inner and outer cylinders, and i?;^ and i?out are 
their angular velocities. 

^From the azimuthal component of the induction equa- 
tion one finds 
B 



B. 







R 



(3) 



with B — fSBoRi^. This parameter /3 therefore denotes the 
ratio of the toroidal field to the constant axial field B^ — Bq. 
This toroidal field B^, is maintained by an electric current 
running along the central axis of strength 



J = 5(3BoRi, 



(4) 



with J in Ampere, Bq in Gauss and Rm in cm. 

We are interested in the stability of the basic state Q 
against axisymmetric and nonaxisymmetric perturbations. 
The perturbed state of the flow is described by the quantities 



Rf] 



B'ji, Bj, + B'j„ Bq + B' 



(5) 



The solutions of the linearized MHD equations are consid- 
ered in their modal representation F' = i^'(i?)exp(i(fcz + 
m(f) + cut)) where F' is any of u' and B' . The dimensionless 
numbers of the problem are the magnetic Prandtl number Pm, 
the Hartmann number Ha and the Reynolds number Re, 



Pm 



Ha 



BqRq 



Re 



(6) 



(RinD)^^^ withZ? 



where fiQ is the permeability, and Rq 

Rout Rin- 

In order to understand the results presented here, it is use- 
ful also to consider the induction equation for the fluctuating 
toroidal field component, 

dB' 



dt 



where ttpoi and .Bpoi denote the poloidal components of u 
and B. Chandrasekhar (1961) cancelled the last two terms 
on the right of this equation, and did not obtain the MRI. It 
occurs if only the last term in is cancelled; but one can 
show that for small Pm it only appears for Rm = 0(10). If 
also the last term in (0 is retained then the resulting Reynolds 
number loses its strong Re cx Pm^^ dependence, and is re- 
duced by 3 orders of magnitude (HoUerbach & Riidiger 2005, 
see also Table|4}. 

3. The equations 

The equations of the problem are given in a similar form as by 
Riidiger & Shalybkov (2004). Lengths and wave numbers are 
normalized by Rq, velocities by rj/ Rq, frequencies by J?in, 
and magnetic fields by Bq. One then obtains 

dP' .m^ ^ 

1- 1—X2 + ikXs 

dR R ^ 



+iRe(cL; + mO)u'ji — 2nRcu'^ -ms? kB'i^ 



-iHa- 



dR r + i?2 



i?2 



0, 



(8) 



M0 — iRc(cj + mfl)u'^ + 2i 



, mB 



-2Rcau'„ + msf—B'^ + iHa^fcSi - i—P' = 0, (9) 
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-ikP' + iHa^^B' + iHa^fcS' = 0, (10) 



du'r, u't, .m . ., , 



(11) 
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(14) 



with P' as the pressure and 

du'^ du' 



X- 



dR 



(15) 



i? di? ' "-^ di?' ^ i? 

The boundary conditions for the flow are u'^ = = tt^ = 
for i? = i?in and R = i?out ■ The boundary conditions for the 
magnetic field are also straightforward; taking the inner and 
outer cylinders to be perfectly conducting, they are simply 



B\ 



R 



0, Xi 



for i? = Rin and R 



consider only the particular radius ratio i) 



= Rout- We wiU 
0.5. 



4. Axisymmetric modes 

Figure |2] shows the critical Reynolds numbers at and beyond 
the Rayleigh line for experiments without toroidal magnetic 
fields, the classical design. We note how Rcc jumps abruptly 
from 10^ to 10® (for liquid sodium, see Riidiger, Schultz 
& Shalybkov 2003). With a toroidal field, however, the so- 
lutions are very different, as shown in Fig. |3l The toroidal 
field strongly modifies the extremely steep line obtained for 
/3 = (Fig. Its inclination is reduced and it starts at lower 
Reynolds numbers. For /3 > 2 we find critical Reynolds num- 
bers of order 10"^. This is a dramatic reduction of the values 
of order 10® which are characteristic for (3 = 0. 

The corresponding Hartmann numbers for these instabi- 
lities are also strongly reduced by the inclusion of a toroidal 
field (Fig.gJ. TablesElandOlshow results for fl = 0.27, (3 = 
to 4, and Pm = lO"'^ (sodium. Table |3 and Pm = 10"® 
(gallium. Table |3}. The first five columns show the nondi- 
mensional quantities indicated. The last four columns show 
dimensional quantities, taking i?in = 5 cm and i?out = 10 
cm. Specifically, Re has been converted to /in, the rotation 
frequency, in Hz, of the inner cylinder. The Hartmann num- 
ber Ha has been converted first to Bq, the strength, in Gauss, 
of the axial field. Next, the product (3 Ha has been converted 
to the strength of the azimuthal field at i?i„. Finally, J denotes 
the axial current, in Ampere, that is required to maintain this 
toroidal field. 

For /3 > 2 we see then that the critical Reynolds numbers 
no longer scale as Pm~^, as they do for (3 = 0, but rather be- 
come independent of Pm. This is the same result previously 
noted by Hollerbach & Riidiger (2005) for insulating bound- 
aries. More quantitatively, we note that for /? > 2 the required 
rotation rates of the inner cylinder are reduced to less than 1 
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Fig. 2. On the Rayleigh line (here fl = 0.25, dotted line) the 
critical Reynolds number scales as Pm""-^, while beyond the 
Rayleigh line it scales as Pm"^, so that for Pm = 10"^ the 
critical Reynolds number jumps from 10'' to 10®. 
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Fig. 3. The critical Reynolds numbers beyond the Rayleigh 
line, for axisymmetric modes (m = 0) with toroidal fields of 
the same order of magnitude as the axial field. Pm = 10"''. 

Hz, the axial fields to a few tens of Gauss, and the axial cur- 
rents to 0(1000) A. All of these values are easily achievable 
in the laboratory. We conclude therefore that this new design, 
incorporating an axial current, is the most promising design 
for obtaining the MRl in a laboratory experiment. 

It is important also to consider the wave numbers k; if the 
vertical cell size is too large the experiment will still run into 
difficulties. We note in Fig.|5]that for non-zero /? fc is reduced 
somewhat, corresponding to a greater extent in z (the vertical 
cell size is given by Sz w ir/k). The increase in cell size 
is not very great though, and therefore should not cause any 
problems. 

Another interesting point in comparing k for (3 = ver- 
sus (3 > is the Pm-dependence, or rather the lack thereof, 
since k is independent of Pm in both cases. The details of how 
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The Hartmann numbers of the solutions given in Fig. 5. The wave numbers fc of the solutions given in Fig. |3l 



Table 2. Characteristic values for i) = 0.5, fi = 0.27, and 
Pm = 10^^. In converting from nondimensional to dimen- 
sional quantities i?in = 5 cm and i?out = 10 cm were used, 
and the material properties of sodium. The magnetic fields 
are measured in Gauss, denotes the toroidal field at the 
inner cylinder, the current J is measured in Ampere. 5R((jj) is 
the real part of the Fourier frequency. 
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1.6 


0.13 


0.05 
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0.15 
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Table 3. 

10-6). 



The same as in Table |2l but for gallium (Pm = 



this comes about are subtly different though. For the classical 
MRI with (3 — the critical wave number is proportional to 
Aii/Va, where Va = Bq/ ^^^p is the Alfven velocity. The 
nondimensional wave number k then becomes 



k oc 



Ha 



(16) 



which does not depend on the magnetic Prandtl number, since 
Re cx Pm^^ and Ha cx Pm^^/^. In contrast, for /3 > 
k is still independent of Pm, but so are Re and Ha, so that 
( fT6t cannot be the relevant balance in this case. For [3 > Q 



the actual wave numbers are much greater than M6\ would 
predict. 

Finally, Fig.|6lshows the frequency of these modes. 
For /3 = this is zero, since the modes are stationary in 
that case. For non-zero /3 stationary modes no longer exist; 
as noted also by Hollerbach & Riidiger (2005), including a 
toroidal field changes the symmetries of the problem in such 
a way that ±z are no longer equivalent, which inevitably 
means that the modes will drift one way or the other in z, 
that is, they will be oscillatory rather than stationary. We see 
though that the real parts of oj are still rather small, with the 
mode-oscillation time exceeding the rotation time of the inner 
cylinder by more than a factor of 10. 
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Fig. 6. The real parts of the eigenfrequency uj of the solutions 
given in Fig.|3] 



5. Nonaxisymmetric modes 



Hollerbach & Riidiger considered only axisymmetric modes 
m = 0. These are usually the preferred modes, both in 
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Table 4. Insulating cylinders: Fields and axial current for 
i] = 0.27 and (3 — 4 (Rm = 5 cm, i?out = 10 cm). The mag- 
netic fields are measured in Gauss, 3^, denotes the toroidal 
field at the inner cylinder, the current J is measured in Am- 
pere. 
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sodium 


1521 


0.07 


16.3 


26 


106 


2657 


gallium 


1521 


0.03 


16.3 


72 


287 


7170 



nonmagnetic Taylor-Couette flow, as well as in the classical 
/3 = MRI. Nevertheless, for a complete analysis the nonax- 
isymmetric modes with m > must also be considered. Fig- 
ure0shows these results, and demonstrates that here too the 
critical Reynolds numbers for the onset of nonaxisymmetric 
modes are greater than for the onset of axisymmetric modes. 
Indeed, including a toroidal field reduces Rcc by far less for 
the nonaxisymmetric than for the axisymmetric modes. For 
/3 > 1 the axisymmetric modes thus occur several orders of 
magnitude before the nonaxisymmetric ones do. These non- 
axisymmetric modes are therefore not interesting from the 
point of view of doing laboratory experiments. 

6. Conclusion 

The central conclusion of this work is as before in Holler- 
bach & Riidiger (2005), that imposing both axial and az- 
imuthal magnetic fields together dramatically reduces the 
critical Reynolds numbers required to obtain the magnetoro- 
tational instability. However, for the insulating boundary con- 
ditions considered there, the axial currents required to gener- 
ate this new toroidal field were at the upper range of what 
could be achieved in the lab. Here we therefore considered 
conducting cylinders, and found that the required currents 
(and Reynolds numbers) are smaller For comparison, Table|4] 
presents the results of Hollerbach & Riidiger for /3 = 4; 
we see that Re = 1521 and Ha — 16.3, compared with 
Re = 842 and Ha = 9.5 here (the last rows of Tables Q 
andO. Switching from insulating to conducting boundaries 
thus reduces both Re and Ha by almost a factor of 2, which 
would certainly help in achieving these toroidal fields in the 
lab. We conclude therefore that implementing this experiment 
with conducting boundaries is the most promising design for 
exploring the magnetorotational instability in the laboratory. 
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Fig. 7. Critical Reynolds numbers (top), Hartmann numbers 
(middle) and wave numbers (bottom) for nonaxisymmetric 
modes with m = 1. 
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